Fine-structure constant variability: surprises for 
laboratory atomic spectroscopy and cosmological evolution of quasar spectra 
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Calculation of the Dirac hydrogen atom spectrum in the framework of dynamical fine structure 
constant (a) variability discloses a small departure in the laboratory from Sommerfeld's formula 
for the fine structure shifts, possibly measurable today. And for a distant object in the universe, 
the wavelength shift of a spectral line specifically ascribable to cosmological a variation is found to 
depend differently on the quantum numbers than in the conventional view. This last result clashes 
with the conventional wisdom that an atom's spectrum can change with cosmological time only 
through evolution of the a parameter in the energy eigenvalue formula, and thus impacts on the 
Webb group's analysis of fine structure intervals in quasar absorption lines (which has been claimed 
to disclose cosmological a evolution). In particular, analyzing together a mix of quasar absorption 
lines from different fine structure multiplets can bias estimates of cosmological a variability. 
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Early theoretical speculations about temporal variabil- 
ity of the fine-structure constant a [ij, |2|, |3j (an excellent 
review is provided by Uzan Q) have assumed urgency 
in view of the claimed cosmological variability in fine- 
structure multiplets seen in absorption in distant quasar 
spectra In the Webb group's analysis it is assumed 

that the spectrum of an atom in an intergalactic cloud 
at redshift z can be gotten from its laboratory spectrum 
by merely replacing the laboratory value of a by that 
at redshift z. This assumption is crucial to the group's 
strategy for reducing experimental errors by analyzing 
simultaneously lines from different fine-structure multi- 
plets of different elements. 

The mentioned assumption takes it for granted that in 
a world where a varies, the form of the Dirac hamilto- 
nian (the most direct tool for determining the fine struc- 
ture splittings) relevant for an electron in an atom is the 
same as in standard physics. Is such an "electromagnetic 
equivalence principle" correct? 

In one paradigm cosmological a variability ,J\ occurs in 
discrete jumps, each precipitated by a phase transition in 
which spontaneous symmetry breaking changes a certain 
vacuum expectation value, and thereby the value of a. 
In this point of view it seems reasonable for the Dirac 
Hamiltonian to retain its textbook form, at least during 
each cosmological interval of constant a. 

Things are different in the much earlier paradigm, first 
contemplated by Jordan 1\ and Dicke where a is 
a dynamical field. Both workers emphasized that in a 
covariant theory, temporal dynamical a variability must 
be accompanied by spatial variability. Concerns, first ar- 
ticulated by Dicke [jj, that spatial variability necessarily 
entails observable violations of the weak equivalence prin- 



ciple have proved groundless, at least within the general 
framework of dynamical a variability proposed by us : 
its violations of the equivalence principle lie beyond fore- 
seeable experimental sensitivities 

As we show here for the hydrogen atom, dynamical 
a variability according to the framework modifies the 
form of the Dirac hamiltonian. One consequence is a 
small deviation of the hydrogen atom's spectrum in to- 
day's laboratory from that given by Sommerfeld's for- 
mula [Eq. (|15|l below]. Dynamical a variability also 
implies nontrivial evolution of atomic spectra with cos- 
mological epoch (time). Estimation of a's cosmological 
evolution from quasar fine-structure multiplets of mixed 
provenance must take cognizance of this effect. 

In its most convenient form, the framework for a vari- 
ability is a theory of constant charges where the permit- 
tivity and the reciprocal permeability of the vacuum both 
equal e -2 ^, where ip is a scalar field [UEl- One conse- 
quence is that the coupling "constant" a of Coulomb's 
law, which is inversely proportional to the vacuum's per- 
mittivity, varies as e 2 ^. Another is that the electromag- 
netic field tensor ft™ an( j its dual obey the equations 
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with j M the electric current. The coupling between e -2 ^ 
and f^ v induces an electromagnetic contribution to the 
source term of the scalar field equation, 
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where the density a is associated with particles and n is 
the parameter of the theory with dimensions of recipro- 
cal energy times charge,. Whenever j M can be described 
as a sum of contributions from discrete particles, a is ob- 
tained by replacing in the product of each particle's 



2 



charge and the corresponding 4-velocity by the deriva- 
tive dmc 2 1 'dip of the particle's mass to. So in this theory 
particle mass is, in principle, a function of the ambient 
tp field. 

In flat spacetimc the physical static solution of 
Eqs. ©-© for a pointlike charge eo in terms of its formal 
Coulomb potential $ = eo/r and its true electric field, 
E^c{f\f 2 J™},i S i 



E 



-e 2 ^V$ = -V(K- 1 tan(K$ + x)) (4) 
sec(K$ + x), (5) 



with x a constant of integration. A second such constant 
has been set to unity by redefinition of the charge 0. 
This restricts our description to epochs for which > 1. 
Now x is evidently equivalent to exp(^ CO sm), the cosmo- 
logical value of at the relevant epoch, which plays the 
role of asymptotic value of e^( r \ Solution ©-© thus 
depends on an adiabatically varying parameter. We as- 
sume that |x| <C 1 throughout the epochs of interest. 

Solution ©-© applies as well to multiple pointlike 
particles at rest 9] with $ replaced by the Coulomb po- 
tential coming from all the charges. According to Ref. 0, 
such a solution is consistent arbitrarily near a source 
only if the latter's mass to is subject to the condition 
k 2 dmc 2 /dtp = —eodip/d<$> with both sides evaluated at 
the relevant particle's position. The theory thus fixes 
each m's dependence on ip. By Eq. © the condition is 



dmc 2 7<9$ = — eo tan 2 («$ + x) 



(6) 



The derivative here gives the change of m reflecting a 
local change of the total Coulomb potential $ caused by 
charges elsewhere (the charge's self-potential, $Sj cannot 
change, of course; since the framework does not allow 
truly point charges, $s < °o @)- 
The integral of Eq. © is 



C-e K 1 tan(^$ + %) + e o® 



(7) 



with C an integration constant. In terms of the free par- 
ticle's rest mass mo (m — > m for $ — > $ s ) at the epoch 
in question, C = mo + f-a^ 1 tan(K$s + x) — eo^s- Now 
from Ref. |9j it follows that k _1 is no more than an or- 
der of magnitude below Planck's energy divided by the 
elementary charge. Thus, for example, even at distance 
10~ 17 cm from an electron (if such accuracy in position 
has any meaning in light of the Zitterbewegung), eo$ re- 
ceives a contribution 18 orders below Planck energy and 
k|$| and k|$s| are both smaller than 10 -18 . Since we 
have no reasons to expect that x is this tiny at a generic 
epoch, we may regard both k|$| and k|$s| as small com- 
pared to x- Thus for a (pointlike) electron with charge 
eo immersed in the Coulomb potential <i>n of a nucleus 
at a cosmological epoch specified by x, we obtain by ex- 
panding Eq. © to first order in K$ n = k($ — $ s ), and 
neglecting in the tangent's argument as compared 
to x that 
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toq + 5m w toq — eoc tan x ' 



(8) 



Consider now the Dirac hamiltonian H for the electron 
in an electric field E. In an inertial frame 

H = (-ihca ■ V + (too + Sm)c 2 (3 + e UI) (9) 

with / the 4x4 unit matrix and (3 and the triplet a the 
four familiar Dirac matrices 12]. By minimal coupling 
the potential U in which the electron moves must be such 
that E = - VU. From Eq. © we infer that 



U = k 1 tan(K$ + x) — K 1 tan(K$ s + x), 



(10) 



at the electron's position, the added constant ensuring 
that U vanishes asymptotically for any x, just as does $n- 
We include the electron's <i>s in because of nonlinearity 
of the theory [see Eqs. 10-©]. Expanding U to 0($ n ) 
and again dropping k$s in comparison with x, we get 



U « sec x • $n- 



(11) 



With help of Eqs. © and <|TT|> . Eq. © may be rewrit- 
ten as 



H = H Q + 8H 
Hq = (—ihca. ■ V + moC 2 /3 
5H = (I - P)e Q tan 2 X ■ $n 



(12) 

e $ n /) (13) 
(14) 



Now the mass of the nucleus (proton) is also subject to 
the same "correction" as in Eq. © : the proton's charge is 
— eo but it feels the electron's potential — $n- However, 
we do not include this mass "correction" in H just as 
we do not add an extra (— eo)(— $n) on account of the 
nucleus to the potential term in Dirac's equation. We 
likewise drop from H the free proton's mass (analog of 
Too) since it cannot affect transition wavelengths. 

Eq. (|13fl is the textbook Dirac hamiltonian for an elec- 
tron (charge eo) in the external Coulomb potential $ n . 
Let a* = eo 2 /he (to be distinguished from the coupling 
"constant" a) . In terms of the principal quantum number 
n = 1, 2, ■ ■ • and the total angular momentum quantum 
number j = |, |, • • • , its eigenvalues, to 0(a* 4 ), are [ll| 



E = m c 



1 



2^2 



a* 
2^ 



1 



3 + 



3 

4?? 



(15) 



The term of 0(a* 2 ) is the Bohr energy level; that of 
0(a* 4 ) is Sommerfeld's fine structure correction. 

In Eq. (|12fl 8H represents the perturbation coming 
from cosmological a evolution. We shall now show that it 
modifies the conventional expression for the wavelength 
A of the transition n — ► n — An, j — > j — Aj to 
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This may be contrasted with the conventional one 
parameter formula coming just from Eq. I|15|) . namely 
Eq. I |lti|l with ce* t— > a everywhere. 
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If 5 denotes the unit-normalized eigenspinor of the 
eigenvalue problem HqE — EqE,, then to first order in 
perturbation theory the cosmological evolution shift in 
the energy level is 

SE = Xtan 2 x (17) 

I = J &(I-0)E e $ n d 3 x (18) 

What is tan 2 \ m Eq. ((T7|) ? Since a/a* is the reciprocal 
permittivity of the vacuum, and exp(2?/> cosm ) = sec 2 Xi 
we may write 



tan 2 x = sec 2 x — 1 = (a — a*)/( 



Thus the permittivity is unity when % = 0. The precise 
epoch in question depends on the cosmological model; we 
defer consideration of this issue to a future publication. 
Our choice precludes study of epochs with a < a*. 

To calculate I one first notes that I — (3 is the matrix 
diag(0, 0, 2, 2) so that with H written as the bispinor 



(6\ 



1 = 2 



J rfv e $n d 3 x. 



(20) 



Now in the usual representation of a in terms of the Pauli 
matrices er, HqE = EqE takes the form ^2] 

- ihccr • V?/ + m c 2 £ + e $ n ^ = £ £, (21) 
-ihca ■ V£ - m c 2 ?7 + e <I>n?7 = E r]. (22) 

By requiring, in spherical polar coordinates {r, 9, ip}, that 
S be an eigenfunction of the z-component as well as of 
the square of the total angular momentum, J (eigenvalues 
mh and j(j + l)/t 2 , respectively), one is led to the two 
possible forms 12] 



or a 



G(r)y m 



i+i/2 



-«5(r)J^r i/2 



(23) 



with 3V m (0, <p) a spinorial spherical harmonic. The two 
forms differ in their parity for given j\ the first corre- 
sponds to the case I = j — | in nonrelativistic theory 
of the hydrogen atom, while the second has £ = j + |. 
If the spinorial harmonics are normalized according to 
/3V mt 3V m sin 9 d9d(p = 1, then the normalization of 
S requires J Q (F 2 + f 2 )r 2 dr = 1 and similarly for G and 
g (as will be clear all four functions can be chosen real). 

In terms of the dimensionless variables x = rmoc/h, 
e = Eo/rriQC 2 and k = j + |, the equations for the first 
form of 5 can be cast into the form 12J 



(e - 1 + a x /x)F - [d/dx + (k + l)/x]f = 0, 
(e + 1 - a*/x)f + [d/dx + (k - l)/x]F = 

These admit the series solution 



(24) 
(25) 



u=0 



f = e> 



(26) 



Consistency between the two relations connecting ao and 
bo requires that 



7 = -l+ Vk 2 ^ 



(27) 



(the negative root is unphysical). In order for the sums 
in Eqs. H26JI to terminate at v = n' = 0, 1, • • ■ (otherwise 
F and / would diverge), one must have 



-1/2 



(28) 



^g) Expanding this to C(a* 4 ) gives the Sommerfeld spec- 



trum, Eq. (|15|l . with n = n! + j + 2 . 

Because $ n is spherically symmetric, the normaliza- 
tion of the spinorial harmonics allows us to write 



X - 



f e $n r 2 dr, 



(29) 



where $ n = —eo/r. If 5 is not normalized to unity, then 
in passing to the integration variable x we must write 



1= - 



2moc 2 a* J Q °° f 2 xdx 
J Q °°(F 2 + f 2 )x 2 dx ' 



(30) 



The exact relativistic wavefunctions |ll| are rather 
cumbersome, but are not needed if all we want is to cal- 
culate I to C(a* 4 ). Thus starting with ao = 1, we used 
Mathematica to work out, for a definite n', all the other 
a v and b v to 0{a* n +1 ) by including in the recursion re- 
lations |12| expressions for 7 and e correct to 0(a»™ +1 ). 
The series in F and / [see Eq. iJHJ] turn out to be of the 
form J2i=o Ui{a*x)a* 21 and a J2i-o Vi{ a * x ) a * 2l i where 
Ui and Vi are polynomials of order n' . Up to normaliza- 
tion, Uo(a*x) is the associated Laguerre polynomial fa- 
miliar from nonrelativistic theory of the hydrogen atom. 
The correctness of the so obtained F up to terms with 
z = l and of the i = term in / was checked by verifying 
that the Hellmann-Feynman theorem [12j, 



dE 
<9a» 



1 da 



|ry\ 2 

= m c 



J co (F 2 + f 2 )xdx 
J °°(F 2 + f 2 )x 2 dx 



(31) 



is satisfied with Eq given by Eq. H15(l . 

A change in the integration variable in Eq. I|30|l to 
y = a*x (essentially the radial coordinate in units of the 
Bohr radius) brings out an extra factor a* in the ratio 
of integrals. Since / carries an overall factor of a*, the 
C(a* 4 ) contribution SE can be computed by including 
only the i = terms of / in the numerator of Eq. I|30|) 
and of F in the denominator, and totally discarding / in 
the denominator. The result of exact integrations is 



1 = 



? 4 
mora, 



1 



3 + 
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2n 



(32) 



The equations for the second form of S in Eq. I|23|) are 
obtained from Eqs. I|26|l - (|28|l by the replacements F 1— > 
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G, f i— > g and k i— > — fe. The coefficients a„ and 6 y 
are thereby changed, but 7 and e are not, except that 
now n' — 1, 2, ■ • • (n is still n' + j + |). Repeating the 
procedure outlined above we find the analog of X to have 
exactly the same form as Eq. (|32|) . One consequence is 
that, at least to C(a* 4 ), the energies of levels with like 
j but opposite parity are still degenerate, as in the exact 
Dirac hydrogen spectrum. 

From Eqs. iJTSJ, (|T7|), (HHJ) and we deduce formula 
dlfijl for the transition wavelength. It shows an uncon- 
ventional dependence on the cosmological "fossil" a* , the 
value of a when the vacuum's permittivity was unity. In 
principle fitting a generic collection of laboratory mea- 
sured transition wavelengths with formula l|16(l should 
allow separate determination of a and a — a* . A reliable 
nonzero result for a — a* (which if the Webb group is 
right, might be of order 10 _5 o;) would establish dynami- 
cal a variability according to the theories represented by 
the framework as opposed to phase transition variability. 
Possibly evidence for this phenomenon lurks in the resid- 
uals of extant measurements of the hydrogen spectrum. 

Turning to quasar spectra we remark that observed 
quasar absorption line wavelengths must first be be cor- 
rected by dividing them by the absorber's redshift fac- 
tor 1 + z. In practice 1 + z is taken as the ratio of 
the observed wavelength of a Bohr transition (n ^ 0) to 
the corresponding laboratory wavelength. The procedure 
precludes separate determination of the a* in the prefac- 
tor of Qlti[l: any difference between it and a is absorbed 
in the mentioned correction, and a* and a can only be 
distinguished by the terms inside the square bracket in 
Eq. (|16(l . The correction also makes it unnecessary to 
allow for the cosmological evolution of mo; its variation 
gets absorbed in the mentioned correction. 

By substracting the expression given by Eq. I|16|) f° r 
A; -1 for the laboratory line from X z , its analog for the 
corresponding redshift-corrected quasar line, we get 



_L 1_ _ m cai 3 (a z - aj) 

X z A; 2irh 



A 



1 



n 3 (j 



1 

2n^ 



, (33) 



where we have replaced a* in the prefactor by ai in ac- 
cordance with our earlier remarks, as well as dropped 
terms of 0((ai — a*) 2 ) and 0((a z — a;) 2 ) in light of the 
assumption |x| <C 1. Contrast Eq. (|33|l with the conven- 
tional formula implicit in the Webb group's analysis, 



1_ _ m cai 3 (a z 



on) 



2nh 



A 



n 3 (j + |) 



3 



(34) 



We observe that formulae l)33|) and 134f) predict exactly 
the same X^ 1 — A; -1 for different members of a single fine 
structure multiplet (constant n and An) because the for- 
mulae differ by a term (oc n~ 4 ) which is constant within 
the multiplet. We conclude that comparison of fine struc- 
ture intervals in a quasar with the corresponding ones in 
the laboratory, with each multiplet taken separately (as 
done in older work by Webb's group Q and its predeces- 
sors) cannot distinguish between the a variability from 
phase transitions and dynamical a variability according 
to the framework. 

However, Webb's group has recently analyzed jointly 
lines from different multiplets and in different species |6j . 
This procedure is sensitive to the difference between for- 
mulae ()33|l and 1|34[) and their nonhydrogen counterparts. 
Whereas the factors in square brackets in the two formu- 
lae sometimes have similar values, there are also glaring 
discrepancies. Often the conventional formula underesti- 
mates the quantity in question by a factor of 2, and for 
the transition 2p 3 / 2 — > 4si m it underestimates it by a fac- 
tor of 6. For the transition 2 3 / 2 — ► 3pi/2 the conventional 
formula overestimates A z _1 — A; -1 by a factor of 31.5 and 
gets the sign wrong. Clearly if a varies dynamically, a 
reliable many multiplet estimate of the cosmological a 
variation from quasar absorption spectra must use for- 
mula (|33J) . And consideration of it is indispensable for 
discriminating between the two mentioned paradigms of 
a variability. 
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